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Abstract 



In this paper we use the corrections to the usual Newton-Einstein 
secular precessions of the perihelia of the inner planets of the Solar 
System, phenomenologically estimated as solve-for parameters by the 
Russian astronomer E.V. Pitjeva by fitting almost one century of data 
with the EPM2004 cphemcrides, in order to constrain some long-range 
models of modified gravity recently put forth to address the dark en- 
ergy and dark matter problems. The models examined here are the 
four-dimensional ones obtained with the addition of inverse powers 
and logarithm of some curvature invariants, and the multidimensional 
braneworld model by Dvali, Gabadadze and Porrati (DGP). After 
working out the analytical expressions of the secular perihelion pre- 
cessions induced by the corrections to the Newtonian potential of such 
models, we compare them to the estimated corrections to the rates of 
perihelia by taking their ratio for different pairs of planets instead of 
using one perihelion at a time for each planet separately, as done so far 
in literature. As a result, the curvature invariants-based models are 
ruled out, even by re-scaling by a factor 10 the errors in the planetary 
orbital parameters estimated by Pitjeva. Less neat is the situation for 
the DGP model. Only the general relativistic Lense-Thirring effect, 
not included, as the other exotic models considered here, by Pitjeva in 
the dynamical force models used in the estimation process, passes such 
a test. It would be important to repeat the present analysis by using 
corrections to the precessions of perihelia independently estimated by 
other teams of astronomers as well, but, at present, such rates are not 
yet available. 

Keywords: Experimental tests of gravitational theories; Modified theo- 
ries of gravity; Celestial mechanics; Orbit determination and improvement; 
Ephemerides, almanacs, and calendars 
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1 Introduction 



Weak-field limits of theories involving long-range modifications of gravity 
recently put forth to address the issues of dark energy and dark matter 
(Dvali et al., 2000; Allemandi et al., 2005; Navarro and van Acoleyen, 2005, 
2006a,b; Apostolopoulos and Tetradis, 2006; Capozziello et al., 2006; Nojiri 
and Odintsov, 2007) are important because such exotic corrections to the 
Newtonian potential allow, in principle, for tests to be performed on lo- 
cal, astronomical scales, independently of the galactic/cosmological effects 
(Capozziello, 2007) which motivated such alternative theories and that, oth- 
erwise, would represent their only justification. In this paper we will show 
how to obtain phenomenologically tight constraints on the viability of some 
of such modified theories by suitably using the latest observational results 
from Solar System planetary motions (Pitjeva, 2005a, b). 

In Section 2 we will consider the four-dimensional models obtained from 
inverse powers of some curvature invariants (Navarro and van Acoleyen, 
2005). After working out the analytic expression of the secular, i.e. aver- 
aged over one orbital revolution, perihelion precession induced by the New- 
tonian limit of such models, we will compare it to the phenomenologically 
estimated corrections to the usual Newton-Einstein precessions of the peri- 
helia of the inner planets of the Solar System. By taking the ratio of them 
for different pairs of planets we will find that the predicted exotic effects 
are ruled out. In Section 3 and Section 4 we repeat the same procedure 
for the four-dimensional model based on the logarithm of some invariants of 
curvature (Navarro and van Acoleyen, 2006b) and for the multidimensional 
braneworld model by Dvali et al. (2000), respectively, by finding that also 
such models do not pass our test. In Section 5 we apply the same strategy 
to the general relativistic gravitomagnetic field (Lense and Thirring, 1918) 
finding that it is, instead, compatible with the ratio of the perihelion pre- 
cessions for all the considered pairs of planets. Section 6 is devoted to the 
conclusions. 

2 The inverse-power curvature invariants models 

In this Section we will address the long-range modifications of gravity ob- 
tained by including in the action inverse powers of some invariants of curva- 
ture not vanishing in the Schwarzschild solution (Navarro and van Acoleyen, 
2005, 2006a, b). From the correction to the Newtonian potential (Navarro 
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and van Acoleyen, 2005) 



aGM fi i,,o . 

v = rw r ' r « r « « 

' c 

where fc is a positive integer number, it follows a purely radial acceleration 

' C 

The length scale r c depends, among other things, on a parameter ji which 
must assume a certain value in order that the model by Navarro and van 
Acoleyen (2005) is able to reproduce the cosmic acceleration (Carroll et al., 
2005; Mena et al., 2006) without dark energy; it is just such a value of 
ji which makes r c ~ 10 pc (k = 1) for a Sun-like star (Navarro and van 
Acoleyen, 2005). Since (Navarro and van Acoleyen, 2006a) 

fc(l + fc) , . 

(6k + 3)2 4fc 3 fc 1 ' 

and r c ~ 10 pc (k = 1), the condition r <C r c for which the expansion in 
r/r c yielding eq. (1) retains its validity is fully satisfied in the Solar System, 
and eq. (2) can be treated as a small correction to the Newtonian monopole 
with the standard perturbative techniques of celestial mechanics. The Gauss 
equation for the variation of the pericentre u of a test particle acted upon 
by an entirely radial disturbing acceleration A is 



duj Vl — 
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■A cos/, (4) 



dt nae 

where a and e are the semimajor axis and the eccentricity, respectively, of 
the orbit of the test particle, n = y GM / a 3 is the unperturbed Keplerian 
mean motion and / is the true anomaly reckoned from the pericentre. The 
secular precession of the pericentre (uj) can be worked out by evaluating the 
right-hand side of eq. (4) onto the unperturbed Keplerian ellipse 

r = a(l — e cos E), (5) 

where E is the eccentric anomaly, and by performing subsequently an inte- 
gration over one full orbital period To this aim, the following relations are 
useful 

( \ — e cos E\ ,_, 
dt = ( J dE, (6) 
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cos/ = - -. (7) 

1 — e cos hi 

Let us start with the case k = 1; the extra-acceleration becomes 

A fe=1 = -^^r 8 r. (8) 
By proceeding as previously outlined and using the exact result 

/ (cos£-e)(l-ecos£0 8 (f£ = [l 28 + 7e 2 (128 + 160e 2 + 40e 4 + e 6 )l 

Jo 64 

(9) 

it is possible to obtain the exact formula 



a=i = -Sv ,GM « I7 ( l - e2 ) 



(10) 

It is important to note the dependence of (cl>) on a positive power of the 
semimajor axis a: this fact will be crucial in setting our test. 

The predicted extra-precession of eq. (10) can be fruitfully compared 
to the corrections to the usual Newton-Einstein perihelion rates of the in- 
ner planets of the Solar System phenomenologically estimated by Pitjeva 
(2005a), in a least-square sense, as solve- for parameters of a global solu- 
tion in which a huge amount of modern planetary data of all types cov- 
ering about one century were contrasted to the dynamical force models of 
the EPM2004 ephemerides (Pitjeva, 2005b). Such corrections are quoted 
in Table 1. They were determined in a model- independent way, without 
modeling this or that particular model of modified gravity: only known 
Newton-Einstein accelerations 1 were, in fact, modeled so that the estimated 
perihelion extra-rates account, in principle, for all the unmodeled forces 
present in Nature. Since July 2005 (Iorio, 2007a), many other authors so 
far used the extra-precessions of the perihelia of the inner planets of the 
Solar System estimated by Pitjeva (2005a) to put constraints on modified 
models of gravity (Gannouji et al., 2006; Iorio, 2006c; Sanders, 2006; Adkins 
and McDonnell, 2007; Ruggiero and Iorio, 2007; Iorio, 2007b), cosmological 
constant (Iorio, 2006a; Sereno and Jetzer, 2006b), various cosmological is- 
sues (Adkins et al., 2007; Fay et al., 2007; Nesseris and Perivolaropoulos, 
2007a,b; Sereno and Jetzer, 2007), dark matter distribution (Iorio, 2006b; 
Khriplovich and Pitjeva, 2006; Sereno and Jetzer, 2006a; Frere et al., 2007; 



1 With the exception of the general relativistic gravitomagnetic interaction, yielding the 
Lense-Thirring effect, and of the Kuiper Belt Objects. 
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Table 1: Semimajor axes a, in AU (1 AU= 1.49597870691 x 10 11 m), and 
phenomenologically estimated corrections to the Newtonian-Einsteinian per- 
ihelion rates, in arcseconds per century (" cy _1 ), of Mercury, the Earth and 
Mars (Pitjeva, 2005a). Also the associated errors are quoted: they are in 
m for a (Pitjeva, 2005b) and in " cy _1 for w (Pitjeva, 2005a). For the 
semimajor axes they are the formal, statistical ones, while for the perihelia 
they are realistic in the sense that they were obtained from comparison of 
many different solutions with different sets of parameters and observations 
(Pitjeva, private communication 2005). The results presented in the text do 
not change if 5a are re-scaled by a factor 10 in order to get more realistic 
uncertainties. 



Planet 


a (AU) 


5a (m) 


w (" cy" 1 ) 


5m (" cy~") 


Mercury 


0.38709893 


0.105 


-0.0036 


0.0050 


Earth 


1.00000011 


0.146 


-0.0002 


0.0004 


Mars 


1.52366231 


0.657 


0.0001 


0.0005 



Khriplovich, 2007), trans-Neptunian bodies (Iorio, 2007c), general relativity 
(Will, 2006; Iorio, 2007a); a common feature of all such analyses is that they 
always used the perihelia separately for each planet, or combined linearly by 
assuming that the exotic effects investigated were included in the estimated 
corrections to the perihelia precessions, and by using their errors to constrain 
the parameters of the extra-forces. About the reliability of the results by 
Pitjeva (2005a), Iorio (2007b) made an independent check by assessing the 
total mass of the Kuiper Belt Objects and getting results compatible with 
other ones obtained with different methods, not based on the dynamics of 
the inner planets. It must be noted that more robustness could be reached 
if and when other teams of astronomers will estimate their own corrections 
to the perihelion precessions. On the other hand, an alternative approach 
would consist in re-fitting the entire data set by including an ad-hoc param- 
eter accounting for just the exotic effect one is interested in. However, such 
a procedure might be not only quite time-consuming because of the need of 
modifying the software's routines by including the extra-accelerations, but it 
would be also model-dependent by, perhaps, introducing the temptation of 
more or less consciously tweaking somehow the data and/or the procedure 
in order to obtain just the outcome one a-priori expects. 

Here we will not use one perihelion at a time for each planet. Indeed, let 
us consider a pair of planets A and B and take the ratio of their estimated 
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extra-rates of perihelia: if eq. (10) is responsible for them, then the quantity 2 



Tab 



co A fa^ 17/2 
a B 



(11) 



must be compatible with zero, within the errors. The figures of Table 1 tell 
us that it is definitely not so: indeed, for A=Mars, B=Mercury we have 

r M aMe = 10 5 ±0.1. (12) 

The situation is slightly better for A=Mars and B=Earth: 

r MaE = 38 ± 3.5. (13) 

An intermediate case occurs for A=Earth and B=Mercury: 

r EMe = 10 3 ± 0.2. (14) 

It is important to note that 

• The uncertainty in Tab has been conservatively estimated as 

5co A 5lo b \ 17 fa A \ 17/2 (5a A 5a B \ 

+ rrRT +-S- TR -A- + -TT (15) 



ST AB < 



u A 



6j B 



|w A | \u B \J 2 \a B J \a A a B , 
by linearly adding the individual terms coming from the propagation 
of the errors in to and a in eq. (4); this is justified by the existing 
correlations among the estimated Keplerian orbital elements 3 

The results presented here do not change if we re-scale by a factor 
10 the formal errors in the semimajor axes (Pitjeva, 2005b) quoted in 
Table 1. The same holds also for the errors in the perihelia rates which, 
however, are not the mere statistical ones but are to be considered as 
realistic, as explained in the caption of Table 1 

The constraints obtained here with eq. (4) are independent of a and 
r c ; should one use eq. (10) for each planet separately to constrain r c , 
it turns out that, for a = 4 x 10" 3 (k = 1), r c < 4.5 AU. Note that 
with such a value the condition r <C r C) with which eq. (1) and, thus, 
eq. (10) were derived, holds for all the inner planets 

For k > 1 the situation is even worse because of the resulting higher 
powers with which the semimajor axis enters the formulas for the 
perihelion rates 



2 It turns out that the multiplicative term depending on the eccentricities has a negli- 
gible effect on our conclusions. 

3 The correlations among the perihelion rates are low, with a maximum of 20% between 
the Earth and Mercury (Pitjeva, private communication, 2005). 
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3 The logarithmic curvature invariants models 



The same approach can be fruitfully used for the model by Navarro and van 
Acoleyen (2006b) based on an action depending on the logarithm of some 
invariants of the curvature in order to obtain a modification of gravity at the 
MOND (Milgrom, 1983) characteristic scale (Sanders and McGaugh, 2002), 
so to address in a unified way the dark energy and dark matter problems; 
in this model the length scale r c amounts to about 0.04 pc for the Sun. The 
correction to the Newtonian potential is 

GMr 3 , . 

Foe— — , (16) 

' c 

which yields the perturbing acceleration 

A oc —^r. (17) 

r c 

By using 

/ (cos£-e)(l - ecos EfdE = -evr(4 + e 2 ), (18) 
Jo 

the secular precession of perihelion induced by eq. (17) is 



M«- -i -(4 + e 2 ); (19) 

' c 

also in this case it depends on a positive power of the semimajor axis; cfr. 
the approximated result by Navarro and van Acoleyen (2006b) for the shift 
per orbit, i.e. 2n (u) jn. 

By taking the ratio of eq. (19) for a pair of planets and comparing it 
to the ratio of the estimated extra-precessions by Pitjeva (2005a) it can be 
obtained 



Aab 



(20) 



The test is not passed. Indeed, for A=Mars and B=Mercury we have 

A M aMe = 30.7 ±0.1; (21) 

the pair A=Earth, B=Mercury yields 

A EMe = 10.6 ± 0.2, (22) 
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while A=Mars, B=Earth A is marginally compatible with zero 

A MaE = 3.4±3.5. (23) 

Note that, even if the real errors in the estimated extra-precessions of peri- 
helia were up to 10 times larger than those quoted by Pitjeva (2005a), the 
pair Mars- Mercury would still be able to rule out the logarithmic model by 
Navarro and van Acoleyen (2006b). 



4 The multidimensional braneworld Dvali-Gabadadze- 
Porrati model 

Another modified model of gravity aimed to explain the cosmic acceleration 
without dark matter is the multidimensional braneworld model DGP (Dvali 
et al., 2000) which predicts, among other things, an extra-rate of perihelion 
independent of the planetary semimajor axis 4 (Lue and Starkman, 2003; 
Iorio, 2005). It is incompatible with the test of the ratio of perihelia as well, 
although less dramatically than the previously examined models. Indeed, 
by defining 



^-1 



(24) 



for A=Mars, B=Mercury we have 

^MaMc = 1.0 ±0.2, (25) 

while A=Earth, B=Mercury yield 

*EMe = 0.9 ± 0.2. (26) 

Errors in the determined extra-rates of perihelion 5 times larger than those 
quoted in Table 1 would allow the DGP model to pass the test. The pair 
A=Mars, B=Earth give a result compatible with zero: 

^MaE = 1.5 ± 3.5; (27) 

the same hold for the other three combinations in which A and B denotes the 
planets with the smaller and larger semimajor axes, respectively. Until now, 
the DGP model was not found in disagreement with the Solar System data 
because the perihelia were used separately for each planet (Iorio, 2006c). 



4 The only dependence on the features of the planetary orbits occurs through a correc- 
tion quadratic in the eccentricity e (Iorio, 2005) which turns out to be negligible in this 
case. 
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5 General relativistic effects: gravitomagnetism 
and the cosmological constant 



It maybe interesting to note that, contrary to the exotic effects induced 
by the modified models of gravity previously examined, the Lense-Thirring 
effect (Lense and Thirring, 1918) induced by the general relativistic gravito- 
magnetic field of the Sun, not modeled by Pitjeva (2005a), does pass our test 
based on the ratio of the perihelia. Indeed, since the Lense-Thirring perihe- 
lion precessions are proportional to a negative power of the semimajor axis, 
i.e. 



must be considered. It turns out that it is compatible with zero for all the 
six combinations which can be constructed with the data of Table 1. This 
result enforces the analysis by Iorio (2007a) in which the extra-rates of the 
perihelia were used one at a time for each planet and linearly combined by 
finding the general relativistic predictions for the Lense-Thirring precessions 
compatible with them. 

6 Conclusions 

In this paper we used the corrections to the Newton-Einstein secular pre- 
cessions of the perihelia of the inner planets of the Solar System, estimated 
by Pitjeva (2005a) in a least-square sense as phenomenological solve-for pa- 
rameters of a global solution in which almost one century of data were fitted 
with the EPM2004 ephemerides, to put tight constraints on several models 
of modified gravity recently proposed to explain dark energy/dark matter 
issues. By using the ratio of the perihelion precessions for different pair of 
planets, instead of taking one perihelion at a time for each planet as done so 
far, we were able to rule out all the considered long-range models of modified 
gravity, in particular the ones based on inverse powers of curvature invari- 
ants by Navarro and van Acoleyen (2005) and on the logarithm of some 
curvature invariants (Navarro and van Acoleyen, 2006b), even by re-scaling 
by a factor 10 the errors in the estimated perihelion extra-rates. The situa- 
tion is less dramatic for the DGP (Dvali et al., 2000) braneworld model since 
if the real errors in the perihelion precessions were, in fact, 5 times larger 




(28) 



the quantity 




(29) 
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than the ones released it would become compatible with the data. Only the 
general relativistic Lense-Thirring effect passed the test. However, it must 
be noted that our results are based only on the extra-rates of perihelia de- 
termined by Pitjeva (2005a): it would be highly desirable to use corrections 
to the secular motion of perihelia estimated by other teams of astronomers 
as well. If and when they will be available our test will become more robust. 
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